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Note on the Integration of a certain System of 
Differential Equations, 

By John Eiesland. 



Professor Craig, in an article which appears in the present volume of this 
Journal and is entitled "Displacements depending on One, Two and Three 
Variables in a Space of Four Dimensions," deduces the following system of 
differential equations : 



■pua + pisy+PziS, 



(i) 



-%f = PnP — PizY + p u h , 

dt 
-£■ = p l3a — p^p _|_ p 3i S f 

-tff = — Pu<* — PiS — PuY » 



which are satisfied by the four groups of direction cosines a,b, c, d; a', b', d, d 1 ; 
a", b", c", d" ; a"', V", d", d>" and admits of the integral 

a 2 + {? + y 2 + & — const. 

We shall first consider the case where the constant differs from zero. By 
dividing the left side by a suitable constant, we may always suppose that 



a 2 + /? 2 + y 2 + 3 2 = 1. 
We now employ the transformation 



(2) 



a = 






. P = 



u 

Vx 2 + 1 



r = 



*/x* + l' 



v x * + i' 
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where x s = X s + (i s + r 8 . The system (1) now takes the form 
-jf=Pu + Pv& —Pisv + % {p u h + p 3i (i + p u v) , 

-^- = Pzi + Pi3% — PkP + v (Pu^ + PuH + Pa*) , 



(3) 



which is a generalization of Ricatti's equation. 

If now we consider a., p, v as coordinates in ordinary space and t as time, 
this system defines an infinitesimal transformation : 



W= \Pu + PnH- — PvP + * CPm* + B4ft + Puv) \ -^ 



+ \PU — Pli*> + i>83^ + M CM< + 1W + Pi?) \ 



dp 

+ \pSi + Pl3^ — 2>23|t* + V (l>14* ■+- P2iP + Psi*) } -^- 



(4) 



(5) 



which is performed on the point % , p , v in the element of time dt . A point 
a,, fz, v is thus transformed from an initial position X , |tt , v into another general 
position X, (i, v by means of a projective transformation which changes with the 
time t. The general integral of (3) will therefore have the form of a finite 
projective transformation 

^ __ «i^o + Vo + Wo + di — <^o + %o + c 3 r + <2 2 

a 4 ^o + *#o + c 4^o + di ' a 4 /L + &4^0 + c 4 v + <Z 4 ' 

v = <h^o + %o + Wo + d 3 

«A) + %0 + <W> + di ' 

where the coefficients a u b it C;, t^, (»= 1, 2, 3, 4) are functions off, and A, , |« , r 
are the coordinates of the initial point and play the role of arbitrary constants 
of integration. It is easily seen that the coefficients a u b t , c { , d t are particular 
solutions of the system (l). In fact, the generai integrals of this are 

a = &!(*! 4- &A + Tcfli + & 4 e?j , 
/? = h^ + h s b s + & 3 Cg + hd z , 
y — h x a z + Tc z b s -f Je 3 c 3 + & 4 d 3 , 
8 = k^i + k^bi + & 3 c 4 + & 4 d 4 , 
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where a lt b lt c x , rf x ; a%, b 2 , c 2 , o\; a 3 , b 3 , c 3 , d 3 ; a it b 4 , o 4 , d 4 form a set of four 
particular solutions of (1), and from this system of integrals we readily pass to 
(5) by taking into account the relations 



a = 



a 



8 ' 



- P _ r 



b ' 



Suppose now that we know a set of three particular solutions of the system 
(3). Let this be a 1( Ag, ^ 3 ; (i lt (i 2 , n 3 ; v lt v z , v 3 , the determinant of which we 
suppose different from zero. We are able to express the general system of 
integrals by means of these solutions. 

An easy calculation will show that this system is 



a* 



a = 



V*5 H- i 



** + 



^•2 



Vxl+l 



i"o + 



V*I + 1 



*0 



A 4 



V^+l 



"0 "f" 



i«0 + 



^0 + 



V^+I V*i + i nu V4 + i u ' vl^TFi 



Pi 



p = 



v«?+ 1 






3<o + 



Pa 



Po + 



Pi 



*\> + 



(li 



Vxl + l ru V^f+i ° ' V^f+T 



P(> + 



V^H-i V*S + i 



v + 



*0 + 



" 2 



("0 + 



"8 



V + 



V«l+ 1 



v*f + 1 u ' V4 + 1 ru ' v«; + i v*f + 1 



V*! + i 



*„ + 



Vxl+ i 



/"o + 



Vx5 + l 



*0 + 



V*j + i J 



(6) 



where 



A 4 = 



Pi Vl 1 
pJ "8 1 
1*3 1>3 1 



^4 : 



A X 


"1 


1 




A 2 


*2 


1 




A3 


"3 


1 





v 4 = 



Aj p, 1 
A a ^ 2 1 
A3 |tf 3 1 



, A 





*1 


Pi 


"1 




Ag 


^1 


"2 




*3 


Ps 


"3 



Returning now to the system (3), we will consider a, fi, v, t as coordinates 

in a four-dimensional space. An integral hypersurface may be represented by 

the equation 

v = pA + o/t + r, (7) 

p, <r, and <r being certain unknown functions of t. Differentiating and substi- 

dv 
~dT 



tuting the values of £± , & , ~ given by (3), we get 
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Pu + Pm*> —PzsP + v (p u h + Pufi + Pup) 

— P \Pli + Pl%P — Plsl* + A (i»14* + PUP + Pw?) \ 

— <t{Pu — Pn*> + P*P + P (Pu% + i>2# + PwP) \ 

„ dp da dt _, 

dt ^ dt dt 

and, since v = pA. + aii + <r is an integral hypersurface, this equation must be 
identically zero for all values of A and [i, when we substitute for v its value 
given by (7). We thus arrive at the following system of differential equations 
satisfied by the functions p, a and r: 

-J- = Pis + i>i 2 <7 + p u r + p (_p 13 p — pttf + ^ 31 r) , 



~ a \ = — Pss—PuP + Pu a + a (-Pisp — P%P + iV*) > 
* = i>34 — J>up— Pu<J + *(Pi3p — Ptfl + iwr), 



(8) 



(ft 

a system analogous to (3). 

Suppose now that we know a particular solution of (8); then the integral 
hypersurface is known, and we shall prove that the system (3) can be reduced to a 
linear system in three variables. Let the surface be written as before 

V = p/l + CjK + T, 

where p, a, t form a particular solution of (8). If we intersect this surface by 
the linear spaces t = const., we get an infinity of planes which stand in a projec- 
tive relation to each other by virtue of the transformation (4). We now intro- 
duce a projective transformation of coordinates, whereby the plane of each space 
t = const, is transformed to infinity. Such a transformation is 

*1 = 7T^ . Pi = k~£ . "i = =r^ • ( 9 ) 

v — pA — ay — t ^ v — pA — ay — t v — pA — ay — * v ' 

Introducing these new variables in (3), taking also account of the equations 
(8), we get the following linear system : 

-fif — —Pl3 — {Pu* + 22>l3p — P i3 <y) *l + (PlZ ~ PvP) Pi + (Pli — Pl3*) Vi , 
--jf = i>23 — {Pli — Pt#) ** — (Pii* — 2p23<7 + i>23p) f-1 + (j?24 + P^) V\ . 

-ft- = —Ps i — (Pu + Psip) *i — (Pu + Pu°) Hi— (p i3 p + Ipu* — Pzsa) v . 



(10) 
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A knowledge of a particular solution of (8) has thus reduced the integration of 
(1) to the integration of the linear system (10), which we shall write in the 
general form 

-ft = «i + «A + aatti + Wi , 

dfii _ 
dt 

dv 1 



dt 



h + W + %i + Vi> 
c i + CjAl + <Wl + C^j . 



J. 



(11) 



We now proceed in the same way as before, writing 

where p x , (T 1( tTj are unknown functions of t. Differentiating and substituting 

the values of — =J- , — ^i , -p- , we get 
o% ctt dt 



Cj + CgA, + c^ + c^ — p x (aj + OgA,! + a 3 ^j + a^) 



— ^l (&i + hK + & 3i"i + Mi) _ ^i 



~- 5 ^i — J 



<& 



Pi 



<ft 






o, 



and this equation must be identically zero for all values of \ and ^ when we 
substitute for v x its value. This gives us the following system of differential 
equations : 

-4r- = c i — («s — c 4 ) pi — btf! — pi (a 4 pi + 6 4 <n) . 






■«: 



'3Pl — (h — C 4 ) <Ti — ffj (« 4 pi + & 4 ffi) » . 



(12) 



dt 

* = Cl — a^ — SjC! — c^j — Tj (a 4 pj + b^^ . 



eft 

Let us suppose that the 6rst two equations, which do not contain the variable t x , 
have been integrated. Substituting the values of o x and 0j in the third it reduces 
to the form 



dt l _ 
lit 



A^ + B, 



which may be integrated by two quadratures. The general integrals p 1? a lt t x 
having been found, all the integral hypersurfaces v x = p^ + a^ + Tj are known, 
and therefore also all the integral curves. We may therefore say : If a particular 
solution of the system (8) is known, the integration of (3) is reduced to the integration 
of a system of generalized Bicatti's equations in two variables and two quadratures. 
33 
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We may now proceed with the system (12), using a method identical with 
the one employed for three variables.* An integral surface may be written 

<ri = $pi + 4. (13) 

where $ and 4 are particular solutions of the equations 

-^-=-a, + (a, — is)* — aW»— ♦ (&<* — &**). 
**= c 2 -c^-(& 3 -c 4 )^-<?)(&^-^). 

Employing the transformation 

p J = h , <r{ = i r , 

we finally arrive at the linear system 



dt 



jgi = A i9 [ + B 2 a[ + 0» 



(14) 



the integration of which can by a method identical with the one employed in the 
case of three variables be reduced to the integration of a Ricatti equation in a 
single variable and two quadratures. We may therefore say, if a particular 
solution p, a, t of the system (8) and also a particular solution $, $ of the system 
(14) are known, the integration of the original system (1) is reduced to the integration 
of a Ricatti equation and four quadratures. 

II. We shall now consider the case where a, (3, y,h is a set of particular 
solutions satisfying the relation 

a 2 + /3 2 + f + b^ = 0, 

and let us suppose that the |)'s are imaginary functions of t. In this case we are 
able to find a transformation by means of which the problem of integration is 
simplified a great deal. 

We put 

a + i{3 _ _ (y + iS) _ x -\ 

y- ih a ~^ ^ ' I (15) 

«— jg _ y + ih \_ \ V ' 

y — iS a + i@ y ' J 



* For a complete discussion of the generalized system of Ricatti's equations in two variables see 
Sophus Lie's " Vorlesungen uber continuierliche Gruppen," p. 778. See also an article by Alf. Gulberg, 
" Zur Theorie der Differentialgleichungen," in Crelle'u Journal, vol. 115. 
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which give us the relations 

_ % l ~ x v — jl l + x v — JL { x + y — _r_ 

x — y 8 ' x — y 8 ' x — y 8 

or, introducing a factor of proportionality p, 

<x=— pi(l— xy), (3 = p(l + xy), y=pi(x + y), 8 = p(x — y). 

Differentiating, we get 

da . ,~ >, dp , . / dx . dy\ 

-w=-^ 1 - x y ) -dt + ^yy-dj + x -dt)' 



dt ~ vt**' dt ^ V 


dt ^ •* dt J' 


%= i( - x + ^t + K 


dx , dy\ 
W + ~dtJ' 


§= c-f)#+p( 


dx dy\ 
~df dt ) ' 


The first two of these equations give us at once 





n . dp da . d8 

— 2i -J- = -= % -4— . 

dt at dt 

Substituting this value of -£- in the last two equations and solving for -=- 

and -Q- , we get 

dt 

. dx _ dy , .dZ_ , (da _ . jg\ 
2l f-dT~W + i ~dT + x \dt % dt)' 

** dt ~ dt * dt ±V \dt dt)' 

Introducing on the right-hand side the variables x, y and p and reducing 
we get the following equations : 

dx ——j(<n — n W I JPu+Pn)i — (Pis+Pzi) 
ST ~ [Pn P3i)X + 2 



+ — (Pi3 + Pu) — ( Pu + Pn) * ^ 



^ = _ iipii - P3i)y+ (PK-P»)i-(P«-P*) 

. — ( PlS — Pto) — (P23 — Pll) * , .2 

f g y ' 

* rf P — « -l (Pu + Pu)i— (Pu + Pzt) _ . ^14 — i>a3 + (j>i3— £m)J „ 
7"5r - " _i,12 + 2 !C+ 2 y ' 



(16) 
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The equations in x and y have the well-known Bicatti's form, and resemble 
those deduced by Darboux in his " Lecons sur la Theorie des Surfaces," vol. I, 
p. 22. We now put 

u = x + — - , v = y + — , 

the above-mentioned equations will take the form 

jfr — — (j>is + p u ) — {Pu + j> 23 ) i 
dt 2 

+ P (Pn — Pad + \ (Pa + P-a) + {Pu + Pn) * f »] ^ . 
<fy* — — (j>is — Ba) — (Pn — -Pn) »' 

* 2 

+ [» (^12 — ^34) + {(#18— JPaO + (P»s — #u)*f «/] /"> 

which may be solved by means of 4 quadratures. We have thus arrived at the 
following result: If a, (3, y, & be a particular solution of (l) satisfying the relation 
a 2 + /3 2 -+- y 2 + b* = , am? */" aZso ^7ie £>'s are imaginary functions of t , then four 
quadratures will be necessary for the complete integration of (l). 

Suppose that the p's are real ; if a; be a particular solution of the equation 

in x, then will T , where x' denotes the conjugate imaginary of x, also be a 

particular solution of the same equation. This is easily seen by changing i into 
— i on the right-hand side and putting x' for x. (Compare Darboux, "Lepons," 
p. 23.) The same holds for the equation in y, so that we get two pairs of par- 
ticular solutions x, — — r ; y, r . But we know that if two solutions of a 

x 1 if 

Ricatti equation are known, only one quadrature is necessary for the complete 
integration of the equation. It follows therefore that when the p's are real 
functions of t , only two quadratures are necessary for the complete integration 
of(l), a, /?, y, 8 being a given particular solution. 



